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1 Introduction 



In the context of 1 + 1-dimensional integrable quantum field theories (IQFTs) , form factors are 
defined as tensor valued functions representing matrix elements of some local operator 0{x) 
located at the origin x = between a multi-particle in-state and the vacuum: 

Ff . . . , 0„) := (0|O(0)|^i, . . . , ^„>J^,...,^„ . (1.1) 

Here |0) represents the vacuum state and \9i, . . . , ^n)™i,...,/i the physical "in" asymptotic states. 
The latter carry indices /ij, which are quantum numbers characterizing the various particle 
species, and depend on the real parameters 9i, which are called rapidities. The form factors 
are defined for all rapidities by analytically continuing from some ordering of the latter; a fixed 
ordering provides a complete basis of states. 

The form factor program for integrable models was pioneered by P. Weisz and M. Karowski 
[U [2] in the late 70s and thereafter developed to a large extent by F. A. Smirnov, who also 
formulated some of the consistency equations for form factors and noticed that they corresponded 
to a particular deformation of the Knizhnik-Zamolodchikov equations [3j. It was found that the 
form factors of local operators can be obtained as the solutions to a set of consistency equations 
which characterize the analyticity and monodromy properties of the form factors and constitute 
what is known as a Riemann-Hilbert problem. It is therefore a priori possible, by solving these 
equations, to compute all n-particle form factors associated to any local field of a massive IQFT. 
This is however a very difficult mathematical problem which has only been fully completed for 
free theories. However, even partial solutions to the problem have proven over the years to 
be extremely useful for the computation of physical quantities. The reason is that all these 
quantities are related in a way or another to correlation functions of local fields. At the same 
time, it is a very well known fact that correlation functions can be expressed as infinite sums 
depending on the form factors of the fields involved. Although such sums can not be performed 
analytically, they happen to converge very quickly in terms of the number of particles of the 
form factors involved. Indeed, convergence is often so good that the knowledge of the two 
particle form factors is enough to get extremely precise results (see e.g. [HE]). This is one of 
the main reasons why, since it first appeared, the bulk form factor program has attracted so 
much attention and has been employed successfully for the computation of correlation functions 
of many integrable models (see |6j for a recent review). 

A natural generalisation of these ideas is their extension to integrable theories with bound- 
aries. The study of IQFTs with boundaries goes back to the works of I. V. Cherednik [7] and 
E. V. Sklyanin [8], where the conditions for a boundary to preserve integrability were estab- 
lished in the form of generalized Yang-Baxter equations which, given the bulk S-matrix, can 
be solved for the reflection amplitudes off the boundary. These pioneering works were followed 
by [HI [Tot flT] ) where the authors were concerned with the computation of reflection amplitudes 
in ATFTs, and by [12) . where the crossing property of the reflection amplitudes as well as an 
explicit realization of the boundary were proposed. The latter work provided in fact the first 
tool for the computation of form factors in the presence of a boundary. The main idea is to 
realize the boundary as a "boundary state" located at the origin of time {t = 0) and expressible 
in terms of the same creation-annihilation operators used to build the bulk Hilbert-space. The 
drawback of this approach is that the boundary state is in fact an infinite sum of particle states 
of the bulk theory, which makes computations very involved, even for free theories |13l I14j . It 
is worth saying that within this approach, only form factors of fields sufficiently far from the 
boundary can be obtained, e.g. form factors of fields sitting at the boundary are not accessible. 
A complementary approach which allows the computation of form factors of fields sitting at 
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the boundary has been proposed recently [15]. It turns out that if one considers the bound- 
ary as a point-hke object sitting at the origin of space (x = 0), then a natural generalisation 
of the bulk form factor program for IQFTs is possible. We will call this the boundary form 
factor program. Within this program the form factors of fields sitting at the boundary can be 
computed for particular models by solving a Riemann-Hilbert problem which is reminiscent of 
the bulk case. Indeed, as for the bulk case, this Riemann-Hilbert problem is very similar to 
the set of equations employed previously for the computation of the form factors of integrable 
spin chains with boundaries |16j-|2Uj. The latter equations result from a generalisation of the 
approach developed by M. Jimbo, T. Miwa and their collaborators [21\ |22| |23] for bulk models. 
The first form factors obtained by means of the boundary form factor program were found in 
[15] for several theories with a single-particle spectrum. An extension of this work followed 
in [23] where higher particle form factors of the sinh-Gordon model were obtained and in |25| 
where a precise counting of the number of solutions to the boundary form factor equations of 
the sinh-Gordon and Yang-Lee models was performed. In this paper we intend to extend this 
program to a class of multi-particle theories: the ^^-affine Toda field theories. 

Affine Toda field theories (ATFTs) have been studied since a long time and have played 
a prominent role in the development of the field of integrable field theories [26^ |27] . Notably, 
the picture proposed by A. B. Zamolodchikov whereby IQFTs may be regarded as conformal 
field theories perturbed by relevant fields was first illustrated on the example of the E'g-ATFT 
[281 129] . From this example and from extended work on classical Toda theory it was expected 
that a different theory should exist for each simple Lie algebra and that this Lie algebraic 
structure would be crucial in the understanding of these models. Subsequently, a lot a work was 
carried out in order to compute the S'-matrices of ATFTs related to each simple Lie algebra. 
Much of this work was reviewed and extended in [30] , where expressions for the S'-matrices of all 
ATFTs can be found. More recently a closed universal formula for the S-matrices of all ATFTs 
which depends solely on Lie algebraic quantities has been obtained |31) . 

Amongst the large class of ATFTs the j4„-case is probably the one for which structures are 
most simple. For example, two-particle scattering amplitudes have poles in the physical sheet 
whose order is at most two (in contrast to other ATFTs where poles of up to order 12 occur). As 
a consequence the pole structure of both bulk |32[[55] and boundary form factors is less involved. 
The 5-matrices of the ^dji-series were first obtained in [34J, where it was suggested that they 
should be given by the S'-matrices of the corresponding minimal ^„-Toda field theories (higher 
rank generalisations of the Ising model, then known as Z„-models) obtained originally in [35] . 
multiplied by a CDD factor which contained the dependence on the coupling constant. 

In the specific context of ATFTs with boundaries, a great deal of work has been carried 
out in the last years to compute the corresponding reflection probabilities and to classify the 
boundary conditions that are consistent with integrability. The first reflection probabilities for 
ATFTs associated to simply-laced Lie algebras were obtained by A. Fring and R. Koberle in a 
series of papers [U [lOl [H] . Notably, the last of these papers dealt for the flrst time with the 
case of dynamic boundaries. Another series of papers by members of the Durham and York 
mathematics groups [6(i\ [57] [5H] achieved the full classification of boundary conditions which 
are compatible with integrability. Later on, many authors have computed further refiection 
amplitudes for various ATFTs (a quite complete list of references can be found in [39]). Given 
the topic of this paper, it is worth emphasizing those contributions which dealt with the ^^-case. 
Besides the work cited above, solutions for all j4„-ATFTs in terms of hyperbolic functions were 
obtained in [361140] . The specific A2-case was studied in detail in [41] . Integral representations of 
the solutions obtained in [36] , as well as solutions for other algebras were given in [12] . Finally, it 
is worth emphasizing that the Ai-case (that is, the sinh-Gordon model) has, independently from 
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other >l„-ATFTs, attracted a lot of interest as its reflection amplitudes exhibit some properties 
that single them out from those of ah other ^„-ATFTs [i3]-[i9]. 

One of the main differences between the solutions in [9l [T0\ [TT] and those in [41^ HOl B2l [36] 
is the fact that whereas the former are invariant under a duality transformation (we will see 
later what this means) the latter are not. However it was found in [39] that these two types 
of solutions can be related to each other in a simple fashion (see section 12. 2p . In fact that 
work provided the first closed solution for the boundary reflection amplitudes which holds for 
all ATFTs and which can be related via ODD factors to all known solutions in the literature so 
far. In this paper we will take the reflection amplitudes computed in [361 132] as starting point. 

The paper is organized as follows: In section [5] we review the expressions and several proper- 
ties of the 5-matrices and reflection amplitudes of the ATFTs. In section [3] we introduce the 
boundary form factor program [15] and compute the minimal one- and two-particle form factors 
for all ^^-ATFTs. We also give a general ansatz for the higher particle form factors, thereby 
characterizing their pole structure. In section [3] we concentrate on the derivation of recursive 
equations for the form factors of A2-ATFT and find solutions for them at the self-dual point, 
B = I, and up to four particles. We show that some of our solutions correspond to the same 
class of fields whose bulk form factors were computed by T. Oota in [32]. In addition, we also 
construct the form factors up to four particles for another class of fields which correspond to 
spin-1 operators of the bulk theory. In section [S] we present our conclusions and outlook. 

2 The ^„-ATFTs : Generalities 

The ATFTs are a family of IQFTs whose exact 5-matrices [33] , -R-matrices [H [Ml dOl (HI 
[^n [^Ol [32] and mass spectrum [35] are completely characterized by an underlying Lie algebraic 
structure related to the simply-laced algebra A„. The theory has n stable particles in the 
spectrum with masses 

avr 

ma = 2m sin— a = l,...,n, (2-1) 

in terms of a fundamental mass scale m and the Coxeter number of An, h = n + 1. If we 
arrange these masses into a vector (mi, . . . ,mn) then this is the Perron- Frobenius eigenvector 
of the Cartan matrix of An- It is common to relate these particles to the nodes of the Dynkin 
diagram of An- In accordance with the symmetry of the latter, the anti-particle of a particle i 
corresponds to the h — i node, as shown in the figure: 



1 2 3 3 2 n = 1 




Figure 1: The Dynkin diagram of the simply-laced algebra An 



2.1 The S'-matrices 



The S'-matrices admit both an integral [501 EI] ^ind a block representation ^34j. The integral 
representation takes the form 



Sab{e,B) = e^v 



dt ^ , . ^ (tQ\ 



(2.2) 
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with 



8smh(f)smh(i(^) sinh ( 



t min(a,fe) 
h 



sinh 



t{h—max(a,b)) 
h 



sinh (^) sinht 



(2.3) 



where B is the effective couphng constant which is a function of the couphng constant /3 which 
appears in the classical Lagrangian of the theory . By construction the S'-matrices are invariant 
under the transformation B ^ 2 — B, which is equivalent to week-strong duality, /? — > 47r//3 of 
the Lagrangian [30]. The point i? = 1 is known as the self-dual point. The kernel above admits 
an alternative representation, based on the identity 



sinh ^ *™™("'^) ^ sinii 



t(h— max{a,fe)) 
h 



Sinh (I) 



E 

p=\a-b\ + l 
step 2 



sinh 



tjh-p) 
h 



(2.4) 



which allows us to bring the S'-matrix into the block product form 

a+b-l 

Sab{9,B)= n W^' 

p = \a — b\ + 1 
step 2 

with blocks 



(2.5) 



or 



{^}6 



{x + l)e{x-l)e 



{x + l-B)eix-l + B)e' 
foo (it sinh (1^) sinh 



exp 



and 



t(2-B) 
2h 



sinhi(^+T) 



sinh i [e 



sinht(l-l) 



t sinh t 



sinh I — 

ITT 



i-KX \ ' 
h ) 



to 



(2.6) 



(2.7) 



From the block form above it is easy to see that these S'-matrices possess simple and double 
poles. As usual in this context, simple poles are related to the presence of bound states in the 
theory which result from fusing processes of the form a + b ^ c whenever a + h + c equals either 



h or 2h. If such fusing occurs the amplitude Sab{(^) will have a pole 9 



ab 



located at 



u. 



ab 




if a + h + c 
if a + h + c 



h 
2h 



(2. 



Closed formulae for the associated S'-matrix residua, which are closely related to the on-shell 
three-point coupling F^^, can be easily obtained. For example, for the process 1 -|- 1 — > 2: 



2 \2 



(Hi) 



-i lim {6 - iul^)Sii{e) 



2sin(f)sm((H^)sin(|^) 



sm 



( {2+B)-K 
\ 2h 



sm 



(4-B)7r 
2h 



(2.9) 



4 



2.2 The reflection amplitudes 



A particular set of solutions for the i?-matrices of ^„-ATFTs was constructed in [36j in a block 
product representation, 

a 

Ra{e,B) = R-a{e,B) = W\\p\\e, (2.10) 

p=i 

with ^ 

II II ^ (^)6> 

The corresponding integral representation was given in [12] 



(2.11) 



' r dt , . ^ (tQ\ 



Raie,B) = R-a{e,B)=exp 
with 

8sinh(i(^)sinh(^(^)sinh(i)sinh(fci 
p,(t,i?)=p,(t,S) = ^ 



(2.12) 



sinh (^) sinht 

83i„h sinh fi»±a^ ^ (^<i^^) . (2.13) 



sinht 

p=i 

A feature of these solutions is that, contrary to the S-matrices (|2.2p . they are not invariant 
under the duality transformation B ^ 2 — B. As a result, they have the property of having 
different classical limits in the weak (B 0) and strong {B 2) coupling limits, namely 

n 

Ra{9,0) = YlSab{0,l), Ra{e,2) = l. (2.14) 

6=1 

This means that in the classical limit the amplitudes Ri{9,0) correspond to a theory with 
fixed Dirichlet boundary conditions whereas Ri{9, 2) correspond to a theory with von Neumann 
boundary conditions |361I371 [38]. How the breaking of duality invariance occurs can be explicitly 
seen by writing the corresponding i?-matrices as 

Raie, B) = Raie, S) J] 5fe J 0, 1 - - J , (2.15) 

where R are the duality invariant i?-matrices constructed in j39j and the identity ()2.15p was also 
found there Q 

*The 7?-matrices constructed in [S^ admit an integral representation of the form p.l2[) with 



8 sinh ( sinh (if) sinh sinh sinh , 



^ ' sinh^ (27^) sinhf 

Notice that this kernel differs from the one given in [39j which contained several typos. 
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3 The boundary form factor consistency equations 

A set of consistency equations whose solutions are the form factors of boundary fields of an 
IQFT was recently proposed in [TS]- There the equations were presented for the special case of 
a theory with a single particle species, but it is trivial to write them for the general case. The 
first three relations 

0\...^J.^^^^^-^^... , n a ^ _ e fa n \T^^\■■■^^i+l^^i■■■( a a ^ ('Q 1 ^ 

[. . . ,tii,tii+i, . . .) — bfj,-^-^^[Pi - UiJ^i)tn [. . . . . .), [6.1) 

j-,0\fJ.i...^J,„_ilJ,„ ,„ „ Q \ _ T) ( Q \T?^\l^l---l^'-n-lt^n/n Q Q \ ( 1 n\ 

-fn (t'l, • • • , t/n-1, — ri^^[tln)rn [t>l, ■ ■ ■ , ^n-l, —f^n) , i-J-^j 

F^^''^''^-''" {01, 62..., en) = R^,{m-ei)FS\^^^^-^-{2^i-ei,92,...^en), (3.3) 

express the fact that exchanging two particles in a form factor amounts to the scattering of those 
particles, which is appropriately taken into account by the two-particle (S-matrix. In addition if 
one of the particles a form factor depends upon is scattered off the boundary with probability 
R^.[Oi), the value of the form factor should remain unchanged, which is expressed by the second 
equation. Finally, the particle might be scattered off the boundary after a 27ri rotation, which 
is expressed by the last equation. For n = 1 and n = 2 the minimal solutions to these equations 
(that is, solutions with no poles in the physical sheet) will play a crucial role in the formulation 
of a general ansatz for the n-particles form factors and details on their computation will be 
provided in the coming sections. 

In addition to these three equations we have two more relations which fix the pole structure 
of the form factors: the kinematic residue equation 

Jim Ff]f^i-^"(^o + ivr,eo,^i,...,^n) = ^ ( 1 - H ^^^,(^0 - ^fc)5/.M..(^o + ^fc) ) 



\ k=l 
-,0\iii...ti. 



{01,..., en), (3.4) 



and the bound state residue equation 



Jim eF^lf''^-''- (0 + ^ + ^ - ^ - e, ^i, . . . , ^„) = ^r^^ifl^ ^i, • • • , On), (3.5) 

where particles a and b fuse to produce particle c which implies that the scattering amplitude 
Sab{e) has a pole at 6 = iu'^f^ with residue given by (FJ^^,)^. 

3.1 The minimal one-particle form factors 

As explained in [15], starting with the boundary reflection amplitudes, it is possible to find a 
minimal solution, Wa{0), to the one-particle form factor equation^ 

Wa{e) = Ra{e)Wa{-9) and Wa{e + ill) = Ra{-e)Wa{iTT " 9), (3.6) 

which takes the form 

Wa{e)=ga{e)ga{t7r-e), (3.7) 

where ga{e) satisfies the equations 

gaiO) = Ra{0)ga{-0) and ga{e)=ga{27ri-e). (3.8) 



^From here onwards we will omit the explicit dependence of the reflection amplitudes, 5- matrix and form 
factors in B. 
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These are the same equations as for the two-particle minimal form factor, with the 5-matrix 
replaced by the reflection amplitudes. Therefore they can be solved in a systematic way, namely 



ga{0) = ga{0,B) = exp 



dt 



t sinh t 



p^{t, B) cosht IH 



TT 



(3.9) 



We can write (|3.9p as an infinite product of Euler's gamma functions by employing the identities: 



and 



We obtain 



dt sinh(at) sinh(/3t)e-'^* _ 1 ^ 
t sinh(tit) 2 ^ 



a+/3+7+ 
2u 



-/3+7+M 



2u 



^ ' 2u 



1 



sinhi 



-2kt 



(3.10) 



(3.11) 



fc=0 



ga{0,B)=ga{9,B) = II 



^sinhl ^-M^ 



cosh 1 f - ^?i-l(e)g(^) 



p=i sinh ^ ^ 



i7r(p— 1) 



COShi(0-^1 nl 



(3.12) 



F gi-BWgW2(g) 



where 



i-ir(t-„ + |;-f)r(fc-„ + ft + f) 

The following identity will be useful for later purposes: 



(3.13) 



+ ^vr)u;a(0) = j| -0 — — . 

p=l eB+2p-2y^^)e2p~B\^^) 



(3.14) 



3.2 The two-particle form factors 



The two-particle form factors are the solutions of the two-particle form factor equations 

Fi\ei,e2) = Sab{e)F2"{02,ei) = Rbie2)Ff {61,-92) = Raii^r - ei)F^\2m - 01,02). (3.15) 

It was shown in [15] that solutions to these equations take the general form 

Ff (01,02) = fab{Ol2)fab{Ol2)Wa{0l)Wb{02)Hyi,y2), with = e'' + 6-'' . (3.16) 

Here 0i2 = 0i — 02, G12 = 0i + 02 and fab{Q) is the bulk two-particle minimal form factor, that 
is a solution of the equations 

fabie) = Sab{0)fba{-O) = fabi^T^i - 0), (3.17) 

with no poles in the physical strip Im(^) G [0,7r]. The function $(2/1,1/2) characterizes the 
particular operator and is a symmetric function of the variables yi , 2/2 j which is the constraint 
given in p]5] for models with a single particle spectrum. In particular, the minimal boundary 
two-particle form factors correspond to $(yi,?/2) = 1- Solutions to (13.170 for ^^-ATFTs were 
first found by T. Oota [32] in the form of an infinite product of F-functions. Such representation 
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can be obtained in the usual way: for an 5-matrix which admits an integral representation of 
the type (12. 2p it is easy to show that a solution to (13.170 is given by 



fabiO) = exp 



1 dt 
2 



^ . ^ ^<^ab{t,B) cosh t{ 1 + - 

Q rsmhr \ it 



(3.18) 



As for the one-particle form factors, we can use ()3.10p and ()3.1ip to express the form factor 
above as an infinite product of gamma functions 



fabiO) 



a+b-1 

n 

p = |a - 6| + 1 
step 2 



{-p)eQl+i-Bi(^)Ql-i+Bi(^y 



(3.19) 



where 



{P)e 



sinhi(^ + M^) sinhi(^ + M|m) 



sinh ^ 



+ 



i-n-jp-l+B) 



sinh ^ 



+ 



in{p+l-B) 
h 



(3.20) 



It is easy to check that the function inside the product above equals the function F™™{9) 
introduced by T. Oota in [32], as it should be. A useful identity involving the minimal form 
factors above is 



fab{0)f-ab{0 + iTT) 



a+b-1 

n 

= |a - fe| + 1 
step 2 



(3.21) 



which was also obtained in [32]. In addition, it is easy to show from (13.190 that 

f-aai29 + nT) 



Combining this with ()3.14p we find that 

mm(a,a) ^mh ( 6 + '^^^^^ 



(3.22) 



faa{29 + i7r)wa{9 + i7r)wa{9) 



n 



sinh 



j7r(p— 1) 
h 



Sinh ^-iZ^ 



Also 



fn{O + f)fn{0-f) 
hi{0) 



p=i sinh [9 + 

sinhl(g + f)sinhi(g-f) 
" sinh i(^ + MizB)) sinh i(^-Miz^ 



(3.23) 



(3.24) 



These identities will be very important in order to bring the kinematic and bound state residue 
equations into a simple form. 



3.3 Higher particle form factors 
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For reflection matrices which have no poles at 9 = m/2 (as is our case here) it is natural to 
make the following ansatz for the n-particle form factors 



k=l 



where }{^^^^'"^" are constants which are invariant under any permutation of the indices ^x, . . . , ;U„, 
and Qn^^^ ^"{yi, . . . , yn) are entire functions of yj = 2 cosh 9i with i = 1, . . . , n. The factors in 
the denominator encode the full pole structure of the form factors. More precisely, the kinematic 
poles are encoded in the {yi + yj) factor, whereas the functions Pab{9) account for the bound 
state poles in the theory. The form of these functions was determined in [32] for the bulk case. 
They are exactly the same in the boundary theory and can be written as, 



P,f,{9) = P-,l{9) = W 2 cosh0-cos(^^^^^^^ , (3.26) 



\a~b\ + l,p^h~l 
step 2 

which is equivalent to Oota's expression (up to constant factors). The value p = h — 1 needs 
to be excluded as it would produce extra poles aX, 9 = in. The main difference in the way 
of encoding the bound state pole structure with respect to the bulk theory is the presence 
of products Pab{dab)Pab{Sab) rather than a single factor Pab{(^ab)- This must be so in order to 
guarantee invariance under change of sign of one of the rapidities. Such prescription was already 
used in [T3] for the Yang-Lee theory for which (j3.26p is simply 2 cosh 9 + 1 (corresponding to 
a = h = \ and h = 2>) and the same holds for the A2-ATFT. 

The ansatz ()3.25p satisfies all form factor consistency equations by construction, provided 
that Qn{yi, ■ ■ ■ ,yn) has the extra property 

From this property we deduce that Qn'^^ ^"{yi, ■ ■ ■ ,yn) niust be a symmetric function of all 
variables related to the same particle type. Therefore we expect that the corresponding form 
factors can be expressed in terms of elementary symmetric polynomials on the variables related 
to each particle type. Plugging the ansatz ()3.25p into the kinematic residue equation we obtain 
recursive relations for the polynomials Qn^^^ ^"'{yi, ■ ■ ■ ^yn) of the form 

Qn+T'"'''"i-y, y,yu---,yn) = Pniy, yi, • • • , yn)Qn (yi, • • • , yn), (3.28) 

where y = 2 cosh 9 and Pniy, yi, ■ ■ ■ , yn) is a polynomial on the variables y, yi, . . . , yn- In addi- 
tion, if two particles a, h fuse to produce a particle c, then also a second set of equations of the 
form 

Q,?-l?''"-^"(y+,y-,yi, • • • ,yn) = Wn{y,yi, . . . ,yn)QSf^'''^"(y,yi, • • • ,yn), (3.29) 



must be satisfied, where y± = 2 cosh \ 9 ± —^j and Wn+i{y,yi, ■ ■ ■ ,yn) is a polynomial of its 
variables. We will now proceed to compute the polynomials Pn and Wn for the simplest model 
of this class: the yl2-ATFT. 
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4 Boundary form factors for 742-ATFT 

We will now write down the equations (|3.28p - (j3.29p for the particular case h = 3 and n = 2. We 
have then a two-particle theory with 1 = 2 and the following fusing processes 

2tt 

1 + 1^2 with = — , (4.1) 

3 

2tt 

2 + 2^1 with = — . (4.2) 

3 

Before writing down the recursive equations for this case, it is convenient to introduce some 
notation. We will write F^^_^"' {{y}m', {y'}n) to denote an n + m-particle form factor where 
m particles are of type 1 and n particles are of type 2. We choose to order indices in such a 
way that all particles of type 1 appear first and all particles of type 2 appear at the end. In 
addition, we will group all variables related to particles of type 1 or particles of type 2 into 
sets {y}m = {yi, • • • , Vm} and {y'}n = {y^, . . . , y^} = {Vm+i, Vn+m}, respectively. Similar 
notation will be used for the polynomials Q. Employing this notation, we can rewrite equation 
(I3:28]l as follows 

Qrn+n+2 ^^{y, {y}m; -y, {y'}n) = Pm+n{y, {y}m, {y'}n)Qm+',C(.{y}m; {y'}n), (4.3) 

where 

m n / m-km-py,m-r , ,m-km-p y,m-r\ 

Pn+m[.y,{y\ni,{y In) = a{y) > > — 7 — r^r^ arapakaaThfyc, 



with 



k,p,r=Oa,b,c=0 \ ^- V+ ^+ ^+ V- ^- 



COsh(26') - COS I 



(4.4) 



u)± = -2cosh ( 6I± — j , = -2cosh( (9±y j , (4.6) 
ry^ = -2cosh(e±^^^^^), A± = -2cosh(e±^), (4.7) 
f)± = -2cosh(e±^^^(^), A± = -2cosh(^±^^^ii±^), (4.8) 

and fj and a are elementary symmetric polynomials on the variables {y}m and {y'}n respectively. 
These polynomials can be defined by means of the generating function, 

n n 

ll{x + yi) = ^x"-^ak, (4.9) 

k=l k=0 

where the subscript k is the degree of the polynomial ak which depends on the n variables 
yi, . . . ,yn- Also from the kinematic residue equations it follows that 

=V3/2l(i7r) ■ (4.10) 

In addition, due to the presence of bound states there is an additional set of recursive equations 
of the form (j3.29p . Employing our new notation, we will write these equations as 

<5mtS"(y+'y->{y}m;{y'}n) = w^m+n(y,{y}m,{y'}n)(5m'rn+i^{y}m;y,{y'}n), (4.11) 



10 



where 

W^+n{y, {y}m, {y'}n) = 9{B) ( 



/ttB 

2 - 2 COS ( — 



k,r,s=0 



T> T_ akCTrCTs, (4.12j 



with 



T± = -2 cosh e± 



ot(1 - 



31/4 



sm 



Tr{2 + B] 
6 



(4.13) 



Notice that g{l) = 1 and that Wm+n does only depend on the variables related to the type 1 
particles. This surprising property is due to ()3.24p and to the following identities 



/ii(2gVi(g + f K(g-f ) 

wi{e)Pui2e) 
fi2{e + f)fi2{e-f) 



1 



2 cos 



/22(e) 



1, 



Pl2{0) = 1, 



(4.14) 



(4.15) 



which hold specifically for the h = 3 case. The bound state residue equation provides also the 
following recursive relations between the constants: 



0\m+2,n 
n+m+2 



rTO\m,n+l 



were we have used the identity 



4/n(27rz/3)2 . 2 ^vr(2 + i?) 
/2i(^7r) = sm 



6 

which is valid for h = 3. Equations ()4.10|) - ()4.16p are solved by 

n+3fc/2 



(4.16) 



(4.17) 



TTO\n,n+3k 
-"2n+3k 



H. 



0\n,n+l+3k 
2n+3fc+l 



H. 



0\n,n+2+3k 
2n+3fc+2 




010,0 



n+3A:/2 



O|0,l 



n+(3fc+l)/2 



O\0,l 
' 



(4.18) 
(4.19) 
(4.20) 



for k gZ+U{0}. 

4.1 Solving the form factor recursive equations 

Finding solutions to (I4.3p - ()4.1ip is very involved, even for the lower particle form factors. For 
this reason here we have made a further simplification and consider only the B = 1 case. We 
will start by introducing some useful definitions. We define the order of the form factor F^]^"" 
as the number [F^]^"'] which characterizes the asymptotic behaviour 



lim_F^|":'"({0 + s}^;{e + 4O~e 



(4.21) 
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where hy 6 + s we mean that each rapidity is shifted by the same amount s. If we were deahng 
with bulk form factors, this would be simply the spin of the operator under investigation. The 
order of the polynomial Q^']!^"' can be defined in a similar fashion. 

It will also be useful to introduce some particular combinations of elementary symmetric 
polynomials, which we will denote by i^I™''^] and Zl™'"! and which are defined by the following 
properties: 

Zi"^+^^-+^]{y,{yU--y,{y%) = 0, K^^^+^^'^y^, y^ , {yU) = 0, (4.22) 

namely, they are the kernels of the equations (I4.3p - ()4.1ip for B = 1. As it turns out, the i^l^'"! 
and polynomials admit a simple closed form 

n n 

^[i,n] ^ akCTl-\ = J2 ^k^i~^, (4.23) 

fc=0 k=0 

and each of the terms in the sums has order n. Here we will also need 

Zt^'^] = 0"2 + <5'l'5"20"l + (T20'\ + Cr\cF2 - 2(T2Cr2 + 0"lCri(T2 + CTg. (4.24) 

These polynomials are in fact the same introduced in ^32j and denoted there by I*"'"] . The K- 
polynomials defined by the property (j4.22p become soon rather complex and no obvious pattern 
in terms of elementary symmetric polynomials seems to emerge. For example: 

i^[2.oi = 3_cj2^(j2, (4.25) 

^[4,0] ^ _3^[2.0](3^2_ (3 ^^^)2)2^ ^[3,01^^(9^^ + ^^) 

-3a3((fTi((3 + 02? - 3^2) - 4^3(3 + a2))(3 + 02) + {a\ + 4(3 + ^72)^)^73) 
+aA{-27a\ + (3 - a2)(3 + ^2)^ + a?(3 + CT2)(27 + ^2(6 + as)) (4.27) 

-criCT3(-9 + al + 3c7i(6cri + erg)) + ((3 + 0-2)(3 + 2^2) + 3cri(3cri + 0-3) - (Ti)aA). 

Similarly, we define the ii't'^'"] polynomials as identical to the ones above up to the replacement 
(T — > CT. Like those, they have order n{n — 1). In this case, not all terms involved are of order 
n(n — 1) but only the leading ones. In fact if we select out only the leading terms we get a new 
set of polynomials which will be related to structures occurring in the bulk form factors: 

Ktl = (4.28) 

^hnlk = ^I'^S + f^2-Kbuik, (4.29) 

^£12 = <ylK'i^l-^l<yAKt^l-aA<j,a,{al + Za,a^-^a,)+al{2al-a,). (4.30) 

These "bulk" polynomials are in fact equivalent (up to a sign) to the polynomials i?i[m,n] iii- 
troduced in |32j and can be expressed as determinants of elementary symmetric polynomials in 
the way described there. We have now all the tools to find explicit solutions to the form factor 
equations. 

4.2 Form factors of "spinless" fields 

In this subsection we will compute all form factors up to four particles of a family of fields, which 
we will denote by O2 and which we will try to identify later. In the title of this subsection we 
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have used the word "spinless" to indicate that the order of the form factors associated to these 
fields is zero and therefore they should correspond to spinless fields in the bulk theory. Indeed, 
it will turn out that this family of operators is nothing but the boundary counterpart of the 
fields whose form factors were obtained by T. Oota [32]. The condition — fixes the 

order of each polynomial Q^']!^" as 

[Q^::'^]=m' + {m + n){n-l). (4.31) 

4.2.1 Vacuum expectation values and 1-particle form factors 

Here we will normalize Qq^^^'^ = 1 and choose 



Qf^l°'^ = ^[o,i], (4-32) 



with ^[0,1] an arbitrary constant. 



4.2.2 Two-particle form factors 

The recursive equations for the two-particle form factors are: 

= o, (4.33) 

Q?l''°(y+,y_) = (y2-3)A[o,i]. (4.34) 

Both equations admit many solutions (in fact, the first equation admits infinitely many). We 
will start by selecting those solutions whose order (I4.2ip is minimal. We then obtain 

Q2^^^'^{yi;yi) = o-i + o-i, (4.35) 

Q^^"'°(yi,y2) = A[o,i](a?-3)+^[2,o]^f''°l, (4.36) 

with ^[2,0] constant. 

4.2.3 Three-particle form factors 

There are three different three-particle form factors that can be obtained from the solutions 
above: Q^^'^'^, Q^^'^'^ and Q^^'^'^. The first one is the solution of the equations 

Q?^''\y;-y,yi) = (y' -3)(3 + yyi -y?)yl[o,i], (4.37) 

Q?"''(yi;y+,y-) = (y^ -3)iA[o,i]{{y + yi? -3) + A^2fi]i3-y^ -yj-yyi)). (4.38) 
The most general solution of these equations with the required order is 

Q3'^''\yi;y[,y2) = ^[l,2]^[^'']i^t°''l - <T2(^[2,0] - ^[0,1])^''''^ + A[0,l]('^2aiai - 3i^[0'2l),(4.39) 

which once again involves a new arbitrary constant ^[1,2]- Let us consider now the three-particle 
form factor involving only type 1 particles. The equation we need to solve is: 

Q^'^'^\y+,y-,yi) = (y^ - 3)(y2 - yff, (4.40) 

and the most general solution reads 

Q?'^'''i{yh) = <^iJ^'''^+A[s,o]K^'''\ (4.41) 
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with 

jM = ^^((3 + ^2)^ -3a?) -4(73(3 + ^2), (4.42) 

and ^[3,0] arbitrary constant. Finahy, since the form factor (I4.35P is symmetric in its variables, 
qO2\0,3 g^ygj^ i^y exactly the same function as Q^^'^'^. 

4.2.4 Four-particle form factors 

There are three four-particle form factors which are related to the solutions above by the bound- 
ary and kinematic residue equations: Q'^^^^'^, Q^^'"'^ and Q^^'^'^- The polynomial Q^^'^'^ is 
the solution to 

Qf'''^'\yi,y2,-y;y) = Q^''^'°(yi,y2)-Po+2(y,yi,y2), (4.43) 
Q?"''(y+,y-,yi;y2) = Q^'^''\yi;y2,y)W,+i{y,y,), (4.44) 

which is given by 

Q?"''({y}3; y'l) = (^[0,1] - ^[2,o])kif^3('^?(-9 + ^2) + 12(3 + ^2) - 4aia3)<Ti 

+j[3,o](3 + ^2)^2 _ 3^[3,o](^^ ^ 4^^)] _ ^jQ_^j(3z[2.i]ai(3a2 - (3 + ^2)') 
+3^3(4(3 + + afai) + ^3(^1 (9 - (72) + 4a3)(3((Ti - ai) + ^1(^2 + ai^i))) 
-(^[3,i]i^t''°^ +^[i,2](i^[''°l +^ij[''°l))^['''^ (4.45) 
where ^[3,1] is a new arbitrary constant. Let us now compute 

(5f''°'^(y+,2/_, 2/1,^2) = Qf''"^'^(y;yi,y2)VF2+o(y,yi,y2), (4.46) 
which is solved by, 

g^2l0'4({y'}4) = ^[0,4]i^f°'^^ - ^[2,0F4(-6<7? + (3 + <72)(9 + <72) " 2(71^3 + 2<74)' 

+ ^[0,1] [3(<T?(-9 + <T2) + 12(3 + <T2) - 4aia3)((3ai + a^f - 3(3 + ,72)^) 

+(74(36(7^ + (3 + <72)^(-45 + (72(24 + ^2)) - 3a?(9 + <72(36 + 7^2)) 
+a3(24a? - 3ai(-l + (72)(9 + ^2) + 4(73(a? - 0-2)) 

+(74(a2(-3 - 4<7? + 7<72) + 4^1(73 - 8(9 + (74)))] 
+ ^[1,2] [-3(3 + <72)(<7f (-9 + (72) + 12(3 + (72))(3(7? - (3 + (72)^) 

-<73(3 + (72)(6(7?(-15 + (72) + (7?(-33 + (72)(J3 + 12(3 + (72)(73) 
-4(71(73(3 + (72)(3(3 + (72)(9 + (72) - (7^) + 12(7f(74(-3 + (72) 
+ (74(3 + (72)(-6(7f (3 + 5(72) + (3 + (72)(81 + 42(72 + <7^)) 
-2(73(74(2,7?(9 - (72) + 3(7i(l + .72)(3 + (72) + 2(3 + (72 + 2(7?)(73) 
+,7l(4,7?(9 - (72) + (3 + (72)(-15 + 7(72) + 8(2,71,73 - (74))], (4.47) 

where ^[0,4] is another arbitrary constant. Finally, we need to solve for the polynomial Q^^'^'^. 



It is the solution of the following three equations: 

Q4^^l'''(y+,y_;yi,y2) = (y' - 3)Q^'^'''\y,y,,y2), (4.48) 

Q'^'^''\yi,y2;y+,y-) = {y' - 3)Qs'^'''{y,yi,y2), (4.49) 

<3f'''^'^(y,yi;-y,y2) = Pi+i(y,yi,y2)Q^^'"^'\yi;y2)- (4.50) 
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Since the polynomials Q^^''^'^ and Q^^'^'^ are identical, the problem is reduced to solving the 
last two equations and imposing invariance under the transformation cr a. The result is 

<3?'''''(M2; {y'}2) = 9<72(3a? -af- 2(72) - 3ai(72(-6<Ti + 2aial + (71(72) 

+ afiai + 2(7i)(-3i^P,o] ^ ^2) ^ 3i^[2,o](3^^(2(Ti + (7i) + a^) 

+ a2(3(i^[2,0])2 ^ 3^[2,0]^^ + ^3) + ^2^9 ^ ^^^^3 + (3 + ^2)^^) 

- 3(71(72(1 + 0-2)((7l + (7i) + (72(3(3(7^ - 2(72)(72 - CTlCrf (12 - (72)) 

+ ^2(3(4 - (72)(3 + (72) + al{9 - (7?(3 - 5(72) - (72(6 + 5(72 - 2(Tl(7i))) 
+ ^[2,0] ^[2,2] (^^^^^j ^[2,0] _ 3(2 _ A^.^o]) - <^2(1 - vl[3,o])) 

- Z[2'2]A[3_oj(9_3^2^^2^2-<T2a2), (4.51) 

with ^[2,2] constant. 

4.2.5 Remarks on operator identification 

In the previous subsection, all our form factors have "descended" either from Q^^''^'^ or from 

Qq^^^'^, so that we had no need to fix Q^^^^'^. One obvious option is to choose Q^^^^'^ = 
and, more generally, to consider form factor solutions which are completely symmetric under the 
exchange of the names of the particles. Such form factors should correspond to a particular type 
of operators which we will call "symmetric". We mean by this, that they should be symmetric 
on the fundamental fields, such as ©12 = 4>i + 4>2 ^^'^ ^12 = 4>idx<p2 + 4>2dx4>i- But there are of 
course many fields which are not symmetric, for example those that depend only either on (f>i 
("type-1 fields") or on 02 ("type-2 fields"). It is natural to assume that for any such fields the 
one particle form factors Q*^^^^'^ (for type-1) or Q^^'^'*^ (for type 2) should vanish. Therefore, 
the solutions which we just constructed could correspond to a "type-2 field" O2 if we impose 
the condition that all form factors descending from Q^^'^'^ are vanishing. In order to identify 
the precise operator our solutions correspond to we should analyze the UV behaviour of the 
two-point function. We can however get some idea by just analyzing the order of the form 
factors obtained. As anticipated before 

[F^iZn = 0. (4.52) 
This indicates that the boundary field to which the solutions correspond to is a spinless field 
in the bulk theory. If we assume, as mentioned above, that this field has p^^^^'^ = 0, then 
a natural candidate is the fundamental field 02 powers thereof. This constitutes an 

infinite countable set of fields, in correspondence with the infinite countable set of solutions 
which the form factor equations seem to have. We should now provide some more reasoning to 
justify our last statement: notice that in all our examples each new form factor involves a new 
arbitrary constant. If we assume that this phenomenon will continue for higher particle form 
factors, then the full set of solutions to our equations will depend on a infinite but countable 
number of arbitrary parameters. One way of understanding this is to assume that the solutions 
obtained correspond in fact to a field which is a linear combination of an infinite but countable 
set of "typc-2" fields, in one-to-one correspondence with the constants ^[TO,n] (since the form 
factor of a sum of fields equals the sum of their form factors). 

4.2.6 Prom boundary to bulk form factors 

It is easy to see that from every boundary form factor solution, a bulk form factor solution can 
be obtained by the simple procedure of shifting all rapidities to infinity and selecting out the 
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leading order terms |25j . A way to see this is to notice that when such shift is performed, the 
boundary form factor consistency equations tend to the bulk equations because of the properties 

lim Ri{e) = lim Wi(9) = lim Sij{e) = lim fij{9) = 1. (4.53) 

6— >oo 6—*oo 6—*oo 9— >oo 

These imply that equations ()3.2p and (13.30 become equivalent to each other and equal to the 
standard crossing relation for the bulk form factors. At the same time, the S-matrix product 
depending on the sum of rapidities becomes 1 in (13. 4p rendering it in its bulk form [21 [3]. 

It is therefore interesting to consider the limit described above for the solutions we have just 



obtained 

QSS'k(yi)=^[o.i], (4.54) 

QSS(M2) = ^[o,iF? + ^[2,o]i^S, (4.55) 

Q^uik (yi; y'l) = c^i + <Ti, (4.56) 

QSulk(2/i;2/'l.y2) = (A[i,2]^S-<^2(A[2,0] -^[0,1]))^'''" +^[0,lF2<Tiai, (4.57) 

^Sui'k (Ma) = ^i^2(^i^2 - 4^3) + A[3,o]^£ik' (4-58) 



<buik ({yis; y'l) = ^liMn - ^[2.o])(^?^3 + ala,){a,a2 - Aa^) - ^[3,i]42z[3'il 

-^[i,2]i^H^t'''^ + (^1^2 - 4a3)(A[o,i]<TifT3(a2 + aiai) - A[^^2]'^ia2Z^^'^^), (4.59) 

<bSk ({y'}4) = ^[o,4]i^K - ^[2,0]^4(<7i - 2aia3 + 2a,f 

+^[0,iF4[4(fTifT3 - o-4)(fTi(T3 + 2CT4) + a^ia^ - 3(TifT3 + Ja^) - 4o-2(it| + aja^)] 
+^[1,2] [-cri(T2crl{aia2 - 4(T3) + ala4,{al - 6(Ji(T3) + (4cT2(<Tf - o-3)fT3 - 8ala1,)a4, 
+{7al - 4ai{aia2 - 4a3))al - 8al], (4.60) 

^Sulk (M2; {y'}2) = (71(720-2 (2(T?0-l + Cri(f7? - 3o-2) - 5(7li^£j!j) 

+(72((7i((7i((7l + 2(7i) + (72) + (Ti((7l(7i - 3(72) + (7l ((7l - 3(Jl)(72) 

+ZP.2l((^p,2]4S - -^2)4'^ + ^[3,0]('^2i^£S + ^242)). (4.61) 

We can now compare these solutions to the ones obtained in [32]. It turns out that, up to 
the three-particle form factors, they completely agree with those when setting q = 1 (which 
is equivalent to i? = 1 in our notation) and matching the constants appropriately. There is 
also agreement for the four-particle form factor (j4.59p but there are some important differences 
between our solutions and those of T. Oota for other four-particle form factors, even though 
their orders agree. The precise differences are: 

<bSk({y'}4) - <Oo;t({y'}4) = 2A[o,i]al(aiCT3 - (74) + (^[0,1] " ^[2,0] )^i^4 (^7^ + 2a4),(4.62) 

and 

<£({y}2;{y'}2) - Q?StJ.{yh;{y'}2) = 3zP.21(i _ A[3,o]){k^'''^ - ^?) 

+<7?(72(0-? + CT2)(3 - E:[2'°1) - a2CTl((7iCT2((72 + CT2) + (Tliafaj + (72(72)). (4.63) 

We believe that Oota's solutions for these two cases must be wrong, as the agreement of all 
other form factors with ours strongly suggests that we are dealing with the same class of fields. 
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4.3 Form factors of "spin-1" fields 



In this section we want to find further solutions to the boundary form factors equations and, as 
a by-product, also new solutions to the bulk form factor equations. We will characterize those 
solutions by the choice 

Q^'°''(y) = %i]y, (4.64) 

where B^q i^ is an arbitrary constant. In addition, we would like our form factors to have 

[FSr"] = l, ^ [Qtt''']=m' + {m + n){n-l) + l, (4.65) 

such that the field O2 can be matched to a spin-1 field of the bulk theory. If we choose Q^^'^'" = 
0, we can identify this field as being of "type-2" again (like dx4>2)- obvious solution to these 
constraints is obtained by multiplying all solutions of the previous subsection by the polynomial 
o"! + ai. If we do so we can identify O2 = dx02, but other less trivial solutions also exist, as we 
show below. 

4.3.1 Two-particle form factors 

A set of solutions, which is compatible with ()4.64p is 

Q2'^''\yi,y[) = <Ti(ai + <Ti), (4.66) 

Q?"'°(yi,y2) = i?[o,i]^i(cT?-3) + (i[2,o]+%o]^i)^"'°^, (4.67) 
where ^[2,0] ^[2,0] constants. 

4.3.2 Three-particle form factors 

The three-particle form factors are given by 

Q?^''\yuy'i,y2) = (^[1,2] + (5[i,2] - i?[2,o])^i + qi,2]<^i)^''''^i^f°''i 

-'^2(i[2,o] + 5[2,o](^i + ai))z[i'2l - S[o,i] [3i^I°'2l<Ti + alaii-e + 2(7? - Saa) 
+CT?(-6 + 2al - 3(72) - ^id-l - 3c7i(-3 + aj + a2 + al)], (4.68) 
Q?"'°({y}3) = (^2 + 3) J[3'01 + a,A[3,o]K^'^'\ (4.69) 
Q?'°''({y'}3) = <r,{al - 2(3 + a2))2 + (5[o,3] + qo.si'^Oi^'"''^ (4.70) 

where all variables ^[a,b]) -^[a,b] C![a,b] constants. Notice that, since (j4.66p is not anymore 
a symmetric function, the polynomials Q'?^^^'^ and Q?^'^''^ are now different from each other. 
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4.3.3 Four-particle form factors 

The four-particle form factors are given by 



Q?"''({2/}3,y'l) = (i[3,l] + B[3,l]ai + q3,i]<Ti)K[3.0]^[3,l] 

-i[i,2]((3 + a2)i-3al + (3 + ^2)') + aiiaj - 4(3 + a2))a3)Z^^^^^ 

-(^[1,2] - i?[2,o])(-4i^f''°l<^i + - 2(3 + a2))'a3)z[3'il 

-qi,2]((3 + CT2)(J[-^'°1 - ai(3a? - (3 + a2)')) + aia.iaj - 4(3 + a2))a3)Zl^'^^ 

+5[2,o]('^i + C7i)[3i^l^'°l(4ai +(7i) - J[3'°la?(3 + C72) - Ua^aia^K^^'^^ 

-aiaiasialiS + ^2) - iaia^)] + 5[o,i] [QXl^'Olaf + 3(7?(3 - al){3al - (3 + (72)') 

-3aiai(18c7^ + (3 + C72)^(15 + 4cj2) - 3a?(33 + I6a2 + 2al)) 

+aj{3{3 + (J2)^(-4 + 3cJ2) + 3aj{-2l + 8^2) + aj{297 + 9^2 - 54al - 8al)) 

-3ai{af + aj{-9 + (T2) + 4(3 + (12)^)^3 - 2(cTf + 4(3 + ^2)^ - 2^2(9 + 2a2))al 

+0-1 (T3( -erf 0-2 + (3 + cr2)^(15 + a2) - 18cr?(7 + 3o-2) + crfo-3(5cr? - 4o-2)) 

-d\a^{2a\ + (J?(9 - 19^2) + 11cti(3 + ^2)^ + 8c7?a3 + 4(3 + a2)a3) 
+a?(cTi(-33 + 9a\ - 10(72)(3af - (3 + a2f) 

+C73(7at - 12(3 + a2f + 3c7?(5 + 3a2) - 4c7ia3))], (4.71) 

Q?"''({2/}2,{y'}2) = 9^i(K[°'^]a? - i^[0'2la2 - <T^)9i^[0'2l(2a? - 3^2)^1 
-3<T2(-18 + 6(7? + al)Gx + 3-ft:[°'2l(Ti(i^[°'2l - a2)(y\ 

+ii'[°'2la2crf + (Ji(27 + 27(T2 + 9(ji + 4a^ - 9cr? - 9a2C^? - o\a\) 
+ai(9K[°'2l + 3(3a? + al) - (9 + 6^2 + 4al)a\ + (K^'^l + a?)af ) 
+'7?(-(<T?(-3 + (7?)) + 3(J2(-1 + g\)) - a\{3ax + C7i(3 - <T2 + a\)) 

+<T?C7l(-5<T^ + 6(7? + <T2(-12 + 5(7?))(72 + (7^(-3ai(2 + <T? - a2)(72) (4.72) 
+(7i(-a20-i(6 + 5(72 - 3<TlCri) + K[°'21(3(7i + C7?) - 6-2(3^1 - C7l)c72) 
+ i[3,0] + qo,3])^l(-3 + o\) + (S[o,3] + qo,3]^l)^2] 

-ZP'2][i^[0'2]K[2'0l(i[2,2] + S[2,2]^i) - i[3,o](^[°''l^i(3 - a\) + a2(<7i + ai)i^[2'°l)], 
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and 



Q?^'°''({y}3,2/i) = (i[o,4] +S[0,4]<^l)i^[°''l 

-(i[2,o] + S[l,2pl)<T4(27 - 6(7? + 12^2 + 2ai<T3 + 2(T4)^ 

[3ai(9(j^(15 + ^2) + 2(3 + (T2)-^(27 + ^2(6 + ^2)) 
-9(T?(3 + a2)(33 + ^2(10 + (72))) + 6(-4(3 + ^2)^ + 3<Tf(21 + (T2))(T3 
-12a?(54 + 27<T2 + 6<t| + <t^)<T3 + <ti(270 + 198^2 + 30a^ - 2al + 3a?(39 + a2))ol 
+4(3a? + 2(3 + a2f)al + 3(Ti(36af - a?(585 + (J2(180 + 7a2)))a4 
+3(Ti(3 + a2)(567 + ^2(231 + ^2(25 + a2)))a4 + 4(Ti(-30 + 3a? - 76-2)6-1^4 
+3(234 + 24af + 114^2 + 6ai - 2al - a?(303 + 80^2 + <T2))a3a4 
+(-4a?(12 + 02) + 3ai(168 + 63^2 + 5a^) - 12a?<73 + 24(4 + a2)<73)<7l] 
+^[1,2] [3(-9<7? - 18a?(3 + (72)^(5 + ^2) + (3 + (72)^(15 + ^2) + 3af(54 + 21^2 + al)) 
-3ai(9af + 6(3 + <72)^(7 + 02) - <7?(135 + 48<72 + ai))a3 

- 39a?(3 + ^2) + (3 + (72)^(15 + a2))al + ai{-a\ + 4(3 + (72))(7i 
+ (3af(-3 + 46-2) + 2(3 + ^2)^(36 + 21(72 + d^l) - (7?(135 + 1536-2 + 39a^ + ai))6-4 
+(7i(-27 - 24(72 - 5(7^ + 2(7?(-9 + 2(72))(73(74 - (5(7? + 4(3 + a2))6-ia4 

+ (-54 + (7?(27 - 4(72) - 3(72 + 5(7^ + 13(71(73 - 7(74)(7S 

+ (^[2,0] - S[l,2]) [-3(71(108(74 - 3(7?(3 + (72)(9 + 02^ + (3 + (72)-^(45 + (72(6 + (72))) 

+6(-54(7^ + 2(3 + (72)^ + (7?(3 + (72)(45 + (72(6 + (72)))(73 + 4(7i(-2(7? + 3(5 + a2))alaA 

+(7i(-108.7? + (3 + (72)(-27 + (-18 + (72)(72))<7i - 4(3,7? + (3 + a2f)al 

-ai(3483 + 12a\ + 25656-2 + 603(7^ + 516-^ + 26-^ - 126-?(99 + 30(72 + ^2))o-4 

+(-48af + 46-?(144 + 396-2 + a^l) + 3(-117 - 57(72 - 3(7^ + (7i))6-3(74 

-46-l(-12(7f - 46-?(73 + 3(4 + ,72)^3 + 6-i(99 + 2(72(18 + 6-2) + 2^4))] 

+qi,2] [-9^1(3(7? - 4(72(3 + ^2)^ - (18 + 21^2 + al) + a? (27 + 99^2 + 45(7^ + 5a^)) 

-3(9a? - 4(3 + a2t - af(27 + 48(72 + -ri) + 4(7?(27 + 45^2 + 15a^ + al))a^ 

-3ai(3(7f - (7?(3 + 13^2) + 12(6 + 5(72 + ^2))^! - (^t - 4(7?a2 + 4(3 + a2f)al 

+(7i(3(7f (-15 + 4^2) - (7?(-729 - 63(72 + 39^1 + al))ai 

+(7i(3 + a2)(-702 + a2(-135 + (72(24 + (72)))(74 + (7i(-9(7? + 8(6 + a2))ului 

+ (^^(-42 + 46-2) + (7?(441 + 84(72 - 5(7^) + 3(-117 - 57(72 - 3(7i + (7i))(73a4 

-al(6-?(-51 + 46-2) - 216-?a3 + 12(4 + 6-2)6-3 + 6-4 (342 + 99^2 -al + II6-4))]. (4.73) 
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4.3.4 Prom boundary to bulk form factors 

The bulk counterparts of the solutions above are: 



Qguik(yi'yl) = -^1(^1 + <^i), (4.74) 

^?Suik(yi'y2) = ai(i?[o,i]a? + i3[2,o]^£S)> (4-75) 

QSi'k(yi>yl>y2) = [(^[1,2] - i?[2,o])^i + qi.ai-^il^t^'^ii^fi - a2B[2,o](^i + <^i)^f'''i 

+ %i][((7i + ai)(a^-2CT?a? + 3(T2fT?) + 2CT2fTi(<T2 + ^i<Ti)], (4.76) 

Cui'k(M3) = a2(aicT2-4cT3)+aii[3,o]i^S, (4.77) 

Qs&y'h) = ^-p3 + ^iqo,3]K]::^, (4.78) 

QtSii{yh,yi) = (^[3,iFi + q3,iFi)4u'ri!^''''^ 

+ (i?[2,0] - i?[l,2])((cx? - 2^2)^3 - 442^1)^"''^ 

- C'[i,2](t^2(f^if^2 - 4fT3) + alai + ai{al - 'ia2)a3ai)Z^^'^^ 

- B[2fi]ai{ai + (Ti)((Jio-2 - 4fT3)(c7icr3 + cjad-i) 

+ %i][-2(cJ?-2cJ2)Vi + CJ3(-(Tf(T2+CT^ + fJl(5cJ?-4fT2)cT3)<Tl (4.79) 

+ (cTl(Ti(9cJ2 - 8CJ?)(19CJ?CT2 - llfT^ - 2cjf)cJ3 - 4(2(7? + a2)(Ti)(Tf 

+ iaial{Wa2 - 9a?) + (73(7(7^ + 9a?a2 - 12ai - 4aia3) + 9i^£|!|<Ti)<7?], 

<32ulk({2/'}4) = %4]<7l-fcK--B[i^2]<7l(74(<7^- 2^1(73 + 2^4)^ 

+ (-B[2,0] - -S[i,2])['5"2('5"l<5-2 - 4(73)(73 - 2(7i(7|(T4 + (72(4(7? + 3(72)(73<T4 

- 4(Ti(2(7? - 3(72)<Tio-4 + 4(4(7? - 3<T2)'T3<7l - 8(7i(,T^ + (74),t1] 

+ %i][-2,T2(^^^2_4^g)^2^(3^^^4_ 3^2(^2 ^2(72)<T3 (4.80) 

+ 4,Ti(3(7? - 7a2)&3)a4 - (<Ti(4(7? - 15(72)(72 + 12((7? - 2(72)(73)<t1] 

+ C'[l,2][-('5"? - 2(72)^(73 - (71 ((7? - (72)<T2<5-4 + (T3((72(4(7i - 5(T?(T2 + 3(72) 

- (7l(9(7? - 8(72)(73)(T4 - (7|((Tl(4(T? - (72)(T2 - 3(7(7? - 4(72)(T3 + ll(Tl(T4)], 

QSul'k ({2/}2> {?/'}2) = (72(7? (^["'2] (7? + ,72(<T2 + <7l(7i)) + 4alai{a2 " (7l(7l)(72 

+ (72(7l((7i — (72(7i)(72 + 3(7i(T2(72 + 3(7i(T2((7]^ — (T2)<72(<5"i + C72) 

+ K[°'2](TiiTi((T2 + (72) - 5(72(71(72 ((T?(a2 - (T?) + (T2(T2) 

- (7i(72((7f(7? - ^2^1) + ZP'2U[3,0](i^[''°1^2((7i + (7i) - K^^M) 

+ i^[°'']z['''k(i[3,0] + qo,3] - 1)^? + '^l(qo,3]^2 - i^[''°li?[2,2])). (4.81) 

Equations (j4.74|) - (j4.81|) provide a set of solutions to the bulk form factor equations hitherto 
unknown. 



5 Conclusions and outlook 

In this paper we have initiated the boundary form factor program for theories with many par- 
ticles, concentrating on the example of A„-ATFTs. We have computed all minimal one- and 
two-particle form factors and provided a description of the pole structure of higher particle form 
factors for these theories. We have then specialize our study to the self-dual point, B = 1, and 
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to the 742-theory, which possesses a pair of particles, 1 = 2, that can be regarded as bound states 
resulting from the processes 1 + 1 — > 2 and 2 + 2 — > 1. The presence of bound states implies 
that besides the kinematic residue equations also the bound state residue equations need to be 
satisfied, which makes the computation of form factors a lot more involved. 

For the j42-case we have obtained all form factors, up to four particles, of two families of 
fields which correspond to spinless and spin-1 fields of the bulk theory, respectively. Indeed we 
have shown that the form factors of our first family of solutions reduce to those obtained in j^32j 
for bulk spinless fields in the appropriate limit. In addition, we have obtained all form factors 
up to four particles of another class of fields which correspond to spin-1 bulk fields other than 
simple derivatives of the previous ones. The bulk form factors of these fields were not known up 
to now and have been obtained as a certain limit of our boundary solutions. 

This work provides further support to the statement that the boundary form factor program 
|15j is a very useful tool for the computation of form factors of boundary fields, even for multi- 
particle theories. The fact remains however that the structure of these solutions is a lot more 
involved than for the bulk case and identifying general patterns is still a challenge, even for 
single particle models. 

We would like to finish by emphasizing that the boundary form factor program for IQFTs is 
still at its early stages of development. The study of more models should provide further insight 
into the mathematical structures of the form factors and applications of the program to the 
computation of correlation functions should be further explored. Focusing on the present work, 
a natural follow up would be to look at ATFTs related to other simple Lie algebras, including 
the completion of the study of the j4„-case which we hope to carry out in a near future. 

Acknowledgments: The author is grateful to Benjamin Doyon and Andreas Fring for their 
insightful comments and careful reading of this manuscript. 
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